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I. INTRODUCTION 



II. THE MODEL 



Negatively (X^) and positively (X+) charged excitons, 
also called trions, have been the subject of intense studies 
in the last years, both experimentally and theoretically. 
The stability of charged excitons in bulk semiconductors 
was proven theoretically by LampertH in the late fifties, 
but only recently they have been observed in quantuip. 
well structures: first in CdTe/CdZnIe by Kheng et aln 
and subsequently in GaAs/AlGaAsJja 

After the initial work bw Lampert charged exci- 
tons in bulk semiconductoraJ as well as in an exactly 
two-dimensional (2D) configurationO were systematically 
studied theoretically. These studies revealed that, due 
to the confinement, the 2D charged excitons have bind- 
ing energies which are an order of magnitude larger than 
the charged excitons in the corresponding bulk materials. 
Apart from these two early studies several works were re- 
centbj|j3ublished on charged excitons in a high magnetic 
field|3i3 where one is allowed to use the single particle 
Landau Level approximation, or in the presence of an 
electric fieldrl 

In order to limit the computational time, the previ- 
ous theoretical calculations used approximations and/or 
simplifications in the Hamiltonian, e.g. replacing the true 
Coulomb interaction by an average interaction, or in the 
wave function, i.e. neglecting the correlation among the 
particles in one or more spatial directions. Because the 
binding energy of the trions is very sensitive to the cor- 
relation between the different particles, it would be in- 
teresting to have a full calculation in order to evaluate 
the approximations that have been made. We present 
here a calculation of the ground state energy for the ex- 
citon and the charged exciton based on the stochastic 
variational approach that fully includes the Coulomb in- 
teraction among the particles (for preliminary results us- 
ing this method see Ref. . The use of the stochastic 
method allows us to handle a big number of variational 
parameters in a reasonable time and to systematically 
increase the accuracy of our solution. 

In the present paper we study the and X^ systems 
in a single quantum well (QW) with a finite height of 
the potential barrier. In the first section we present the 
Hamiltonian of the problem. In the second section we 
discuss the dependence of the charged exciton correla- 
tion energy on the well width and on the hole-mass. We 
compare our results with those of Stebe et alSJi, where a 
variational technique with a 66-terms Hylleras trial wave 
function was used. In this section we also present our 
results for the binding energy of the X~ and X+ and 
we discuss the pair correlation functions and the proba- 
bility density of the system. Our results are compared 
with available experimental data from the literature. In 
the last section we summarize our results and give our 
conclusions. 



In this section we present the Hamiltonian describing 
a charged exciton and we discuss the technique which we 
used to solve it. In particular, we focus on the Hamil- 
tonian of a negatively charged exciton. The positively 
charged exciton Hamiltonian can be easily obtained from 
the X~ by replacing the electrons by holes and the hole 
by an electron. 

The Hamiltonian of a negatively charged exciton in 
a quantum well is, in the effective mass approximation, 
given by 



H = Ti,+ T2e +Th + Vc + Vie + V2e + Vh, 



(1) 



where le, 2e indicate the electrons and h the hole; Vie^ 
Vh are the quantum well confinement potentials; Ti is the 
kinetic energy operator for particle i, 



T, = 



Pi 
2m, ' 



(2) 



with TOj the mass of the i-th particle; Vc is the sum of the 
Coulomb electron-electron and electron-hole interactions. 
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Vc = — . _, _, _, . 
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(3) 

with e the elementary charge and e the static dielec- 
tric constant. For a GaAs/A^Gai-jjAs quantum well 
the heights of the square well confinement potentials are 
Vie = 0.57 X (1.155a; + O.STx^) eV for the electrons and 
Vh = 0.43 X (1.155a; -I- 0.37a;2) eV for the hole. 

The Hamiltonianpis then solved using the stochastic 
variational method.E3 The trial function, for the varia- 
tional calculation, is taken as a linear combination of 
"deformed" correlated Gaussian functions, 



0o(^le,^2e,^h) 



K 



^•no(^le, ^2e, ^h) = ^ < CXp 



^ ^ A-nijkOTikT jk 



i,je{le,2e,h} 
kG{x^y^z} 



where r^^ gives the position of the i-th particle in 
the direction k] A is the antisymmetrization operator 
and {Cno, Anijko} are the variational parameters, and 
^(1, 2, 3, ) is the spin function. The '0' in Eqs. (|) refers 
to the ground state. Note that in contrast with the "clas- 
sical" correlated Gaussians here, the parameter Anijko 
which expresses the correlation among the particle i and 
the particle j in the direction fc, is allowed to be different 
from the parameter A^^^f./^ which couples the same two 

particle i and j in a different direction k . This additional 
degree of freedom in the calculation allows us to take into 
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account the asymmetry introduced in the 3D space by the 
presence of the quantum well. The dimension of the ba- 
sis, K, is at first increased until the energy is accurate 
to the second digit, a typical value of K in this calcula- 
tion is 300, and then is refined to increase the accuracy. 
The refinement is made by replacing the n-th state with 
a new state, i.e. with a state built using new parameters 
Cno , Anijko in such a way that it lowers the total en- 
ergy. The process is reiterated multiple times for all the 
K states, until the energy reaches the desired accuracy. 
One can get faster convergence by taking into account 
the cylindrical symmetry choosing Anijxo — j^nijyn- 

III. THEORETICAL RESULTS 

The correlation energy of a charged exciton is defined as 

Ec(X-)-E(X-)-2Ee-E,., (5) 
Ec(X+) = E(X+)-2E,, -Ee, (6) 

with E(X^) the energy level of the charged exciton 
and Ee and E/j the energy levels of the free electron 
and hole, respectively, in the quantum well. Thus, 
Ec is the energy due to the Coulomb interaction be- 
tween the charged particles. We discuss here the re- 
sults obtained for a GaAs/Al^jGai-ajAs quantum well 
with X = 0.3, where the value of the masses used are 
nif, — 0.0667mo, rrihh = 0.34mo, i.e. GaAs masses, with 
Ry — h/m^aB = 5.79 meV, the donor effective Rydberg, 
and aB — h^s/mee^ = 99.7 A, the donor effective Bohr 
radius. The results for the correlation energy of the exci- 
ton and the X^ are shown in Fig. ^ and compared with 
the theoretical results of others. We observe that the cor- 
relation energy of the exciton, Ec(X)=E(X)— 2E?i— Eg, 
(dotted line in Fig. increases in absolute value for 
well widths up to L= 30A, were it reaches a minimum 
Ec(X)= -11.7 meV. For L< 30A and with decreasing 
L the exciton correlation energy decreases in magnitude 
due to the fact that the electron, and to a lesser ex- 
tent, the hole wave functions spill over into the barrier 
material of the quantum well. Consequently the exci- 
ton becomes more extended in the z-direction and the 
Coulomb interaction among the particles is diminished. 
At L=0 we obtain E(7(X)=:— 4.80 meV which compares 
with the correlation energy of an exciton in bulk GaAs, 
i.e. E^^(X)=-4.84 meV. For L>30A and increasing L 
the correlation energy decreases in magnitude with L, 
which is due to the fact that the electron and the hole 
are more extended in the z-direction. In the limit L — > oo 
we recover the 3D exciton in bulk GaAs. 

The correlation energy of the negatively charged exci- 
ton has the same qualitative L-dependence as the ex- 
citon. It reaches a minimum at about 30 A with 
Ec(X-)=-13.2meV. For L>30A it proceeds almost par- 
allel to E(7(X), in the region shown. For the X~ we obtain 
E^^(X~)==— 4.95 meV as 3D correlation energy. 



We compare our results 

for the X in GaAs/Alo.aGao.yAs to the ones reported in 
Ref. ^ (short-dashed curve in Fig. |l]), which are derived 
using the theoriuin Ref. |lj. The theory of Andreani 
and Fasquarelldlil also includes the non isotropy of the 
masses, the non-parabolicity of the conduction band and 
the dielectric constant mismatch. Moreover the values 
for the heights of the potential barrier used are slightly 
different with respect to ours. However, a comparison be- 
tween the two calculations can still be made. We observe 
that our results and the ones in Ref. |l^, are very close 
in the range 70-120 A. For L<70A, our values for the 
correlation energy are much smaller, in absolute value, 
than the ones reported in Ref. |l3|, which is due to the 
band non-parabolicity which is known from Ref. ^ to 
be the major factor responsible for the steep decrease of 
the correlation energy at small quantum well widths. We 
want now to estimate, although in a naive way, the effect 
of band non-parabolicity on our results. From Ref. |lj 
we estimate that for a quantum well L=30A, and a;=0.3 
the parallel mass of the electron is about O.OSmo. If we 
consider, in a very simplified picture, that: a) the con- 
tribution to the correlation energy strongly depends on 
the mass along the growth direction, namely the confine- 
ment energy, which has been subtracted out in the cor- 
relation energy, and that b) the energy of the exciton is 
largely dominated by the mass of the electron, through 
the change in the Rydberg which re-scales the energy. 
Such a procedure gives Ec=—^i-Q meV, for the case of 
a quantum well of width L=30A, which compares very 
well with the value —14 meV given for a quantum well 
of width L=27 A in Ref. |l|. 

Both for the cases of an X and of an X^ in a 
GaAs/Alo.sGao.vAs quantum well we compare our re- 
sults to the one obtained in Ref. Our calculation 
gives qualitatively the same correlation energy both for 
the X and pss: the X~ as compared to the one given by 
Stebe et alx3 However while for the exciton we find that 
the correlation energy is lower than the one obtained in 
Ref. |ll|, thus indicating that the Coulomb correlation is 
more fully included in our approach, for the negatively 
charged exciton our approach gives a higher correlation 
energy (about 4% for L=100A). The latter can be under- 
stood as follows: in Ref. ^ the Coulomb potential along 
the z-direction was approximated by an analytical form 
and Hylleras-type functions were used for the wave func- 
tion. They calculated the Coulomb potential matrix be- 
tween any two basis states (si, s2) by integrating it over 
the /9-plane thus obtaining a potential matrix Vsi^s2{z). 
Then Vsi^s2{z) was replaced by the analytical expression 

--il{6+\zie-Zh\)-l/{6+\zie-Zh\)+l/{l3+\zie-Z2e\) 

where 7, 5 and (3 were determined in such a way that 
it reproduces the correct behavior of the Coulomb po- 
tential matrix in the limit of zero and infinite distance 
between the particles in the z-direction. This approxi- 
mation leads, as the authors of Ref. |ll| noted, to an error 
in the exciton correlation energy which was estimated to 
increase its absolute value by approximately 5%. Our 
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present results for the exciton energy are about 8% lower 
than those of Stebe et aZ.EJ For the charged exciton en- 
ergy the authors of Ref. |ll| did not report an estimate 
of the error which was introduced through the approx- 
imations made. However, we find a smaller correlation 
energy of about 4% as compared to those of Ref. We 
think that this result is not in conflict with the one ob- 
tained for the exciton. Indeed, while in the exciton case 
in Ref. ^ only the attractive interaction between the 
electron and hole was underestimated, for the charged 
exciton case the repulsive interaction between the elec- 
trons will also be underestimated. The difference is that 
the former interaction has the effect of increasing the 
bonding of the particle while the latter has the effect 
to diminish it. Our result indicates that a larger error 
is made in the electron-electron repulsive interaction in 
Ref. |ll| as compared to the error in the electron-hole at- 
tractive interaction. _ 

The approximation by Stebe et aZO consisted in aver- 
aging the wave function in the xy-plane in order to find an 
effective Hamiltonian describing the exciton and the trion 
in the z-direction. This is similar to an adiabatic approx- 
imation which is valid when the motion in the xy-plane is 
faster then the one in the z-direction. We believe that it is 
more natural to do the reverse and average over the parti- 
cle motion in the z-direction which is due to the quantum 
well confinement and which will be much faster. Such 
an approach is equivalent to neglect the particle-particle 
correlation along the z-direction which we expect to be 

valid when i?e(;i) >> E^' . For both the exciton and 
the trion this relation is satisfied for L<150 A. Averaging 
Eq.(|I]) over z we obtained an effective 2D Hamiltonian, 
in which the effective Coulomb potential was replaced by 
the analytical form e^/e[A -I- (pi — Pj)^]^^^, where A was 
obtained by fitting this analytical form to the numerical 
results for the effective Coulomb interaction. The corre- 
lation energy for the exciton and the charged exciton is in 
this case lower than the one we obtain with our more ex- 
act calculation presented above, e.g. in the frame of the 
model presented in this paper we find Ec(X)= — lO.lmeV 
and E(7(X~)= — 10.9meV for a 100 A wide quantum well, 
and Ec(X)= -10.7meV and Ec(X-)= -11.6meV for a 
80 A wide quantum well, while using the screened 2D 
Coulomb potential we found E(7(X)= — 10.4meV and 
Ec(X^)= — 11.4meV for a 100 A wide quantum well 
and Ec(X)= -llmeV and Ec(X-)= -12.2meV for a 80 
A wide quantum well. Consequently, such an approach 
leads to larger correlation energies and also to slightly 
larger binding energies. 

In Fig. |l| we also report the result of a simplified model 
(open diamonds) for the study of the energy of a trion 
which we proposed in Ref. |l^. 31is model is derived from 
the one used for a D~ system,E3 and it assumes that the 
hole is fixed at the center of the well, i.e. it has an infinite 
mass. The effect of the hole is reflected in the renormal- 
ization of the mass of the electron, i.e. mg is replaced 
by the reduced mass /i — rriemh / {nie + rrih). This model 



gives for the correlation energy of the charged exciton 
results that are, at first- surprisingly close to the one ob- 
tained by Stebe et aZ.,0 at least down to well widths of 
about 40 A. This seems to suggest that the procedure of 
averaging the potential in the plane adopted in Ref. ^ 
is almost equivalent to localize the hole in the p-plane. 
For smaller well width the magnitude of the correlation 
energy becomes much larger as compared both to our 
present result and to the one of Ref. |ll]. As shown in 
Fig. the result we obtain in the L=0-limit is dramat- 
ically different from the one found in the present work: 
Ec=— 5.2 meV. The reason is that for small well widths 
the hole may no longer be considered as a " fixed" parti- 
cle and the penetration of the hole in the barrier can no 
longer be neglected. 

To prove further the accuracy of our calculation and 
to check the quality of our wave function for the X^ we 
calculated the virial which is defined as 

{(t)N\w\(j)Ny 

with T the total kinetic jeaergy operator and W — 
Yj^^^TidV/dri. It is knownta that for a system of par- 
ticles interacting through the Coulomb interaction this 
quantity has to be 1 for the exact wave function. We 
obtained a value of 0.999 for almost all the quantum well 
widths studied, which suggests that our wave function is 
well chosen. 

Next we investigate the effect of taking a different mass 
of the particles in the well (GaAs) and in the barrier 
(Al2;Gai_a;As) material, which is expected to be impor- 
tant in the narrow well regime where the electron and 
hole wave functions penetrate into the barrier (see inset 
of Fig. ^). The values for the GaAs-masses, i.e. the 
masses for the electron and the hole, are taken equal to 
the one used in the previous calculation. The values for 
the masses in Ala;Gai_a;As are m*^ = 0.067 -I- 0.083a;, 
m^j, = 0.34 -t- 0.42x, where x indicates the percentage 
of Al present in the alloy. If we assume, as a first ap- 
proximation, that part of the electron and the hole wave 
function is in the quantum well and the rest is in the bar- 
rier we may take the total effective mass of the electron 
and the hole as given by 



rrii rtiiw niib ' 

where , mn, are the masses of the i-th particle in the 
barrier and in the well, and Piw , Pib are the probabilities 
of finding the z-th particle in the well and in the barrier, 
respectively. The results of this calculation are shown in 
Fig. 1^ for X — 0.3. The correlation energy increases in 
absolute value and this is consistent with the fact that 
the effective masses are now larger. The effect of the 
mass mismatch is important only in the narrow quantum 
well regime, i.e. L< 40 A, where it leads to a substan- 
tial increase of the magnitude of the correlation energy. 
In the L=0 limit we obtain now Ec(X^)= — 7.5meV 
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which compares to the 3D correlation energy of a trion 
in Alo.3Gao.7As which we found to be E^^(X")= -6.6 
meV. The minimum of the correlation energy is now ob- 
tained at L=17 A. 

We also studied the dependence of the total energy on 
the hole mass for a 100 A and a 200 A wide quantum 
well. The result, reported in Fig. ^, shows that the total 
energy decreases as the hole mass increases. The energy 
of the negatively charged exciton approaches the energy 
of the D~ in the same quantum well from above and they 
become practically equal when mh/rrig > 16 for the 200 A 
quantum well. Note that for large values of the hole mass 
the X+ energy is practically parallel to the one of the . 
In fact, if the hole mass is large its confinement energy 
contribution to the total energy is negligible, and the 
difference between the X+ and X~ total energies is just 
the confinement energy of one electron, which of course 
does not depend on the hole mass. 

Next we studied the correlation energy of the positively 
charged exciton. In Fig. ^ we plot the correlation energy 
of the X~ and X+ systems as function of the well width. 
Note that the correlation energy of the X+ is equal to 
the one of X~ (within the numerical accuracjy). This is 
in agreement with recent experimental dataEI where the 
binding energy of the X+ was found to be equal to the 
one of the X~ . In fact we have 

Ec(X-) - Ec(X+) = E(X-) - Ee+ E;,- E(X+) 

= [E(X-) -E(X)-Ee] - [E(X+) 
= Eb(X-) - Eb(X+), 

where Es(X^) is the binding energy of a charged exciton 
system referred to the one of the exciton plus one free 
electron (hole) system, 

Es(X±)=E(X) + E^(,) - E(X±), (9) 

where E(X) is the energy of the exciton, Eg(^) is the en- 
ergy of the free electron (hole) and E(X='=) is the charged 
exciton binding energy. Consequently, if the X~ and the 
X+ correlation energy are the same, the corresponding 
binding energies will also be the same. 

Last we study the wave function of the negatively and 
positively charged exciton and the correlation between 
the different particles. The pair correlation function, 
gfj^{r) = {S {r — |T*i — T*j|)), for a 100 A wide quantum 
well is shown in Fig. ^. This function gives the prob- 
ability to find particle i and particle j at a distance r 
from each other. Notice that glf^{r) is the same for both 
X^ and X+ (dashed curve in Fig. ||) and in both cases 
the electron and the hole tend to be close to each other. 
A similar result is obtained for the exciton (dot-dashed 
curve in Fig. ||) . The fact that the intensity of the corre- 
lation function for the exciton is higher than the one of 
the charged exciton is a direct consequence of the normal- 
ization of the wavefunction to one. The situation is very 
different for the correlation between particles having the 
same charge. For the X~ electrons g^^{r) (dotted curve 



in Fig. H) shows that the two electrons avoid each other 
at small distances and have the highest probability of sit- 
ting at a distance of 25 A« 0^/4. For the holes in X+ 
gfi^{r) (solid curve in Fig. ||) shows that the two holes 
avoid each other at small distances and have the high- 
est probability of sitting at a distance of 80 Aw 40^/5, 
thus farther from each other than the electron-couple in 
X~. However the average distance, (|T*i — T*j|), of the 
two electrons in X~ , does not differ much from the aver- 
age distance between the holes in X+. We found 250 
A and 216 A respectively. The average distance be- 
tween the electron and the hole is 150 A and is found 
to be the same in the X^ and in the X+. In the in- 
set of Fig. H we show the 2D-pair correlation function, 
gffip) = {S[p — \7^i ~ 1^ for a 100 A wide quantum 
well, where the same curve conventions are used as for 
the 3D pair correlation functions. These 2D correlation 
functions express more clearly the Coulomb correlation 
between the particles. In the 3D correlation functions 
the z-direction is still involved . In this direction the 
quantum well potential forces the particles towards the 
middle of the well. As a consequence all 2D correla- 
tion functions are more spread out as compared to their 
3D counterpart. The peaks in the electron-electron and 
hole-hole correlation functions are shifted towards larger 
distances. The average distances in the p-plane of the 
electrons in X~ and of the holes in X+ is {\~]f i — ~7^i\) — 
249 A and 214 A, respectively. This result differs only 
— tE(K)fcw E^^stroms as compared to the 3D result, sug- 
gesting th^-jthe charged exciton, for L=100 A, is almost 
bidimensional. 

We now look at the 2D correlation function for differ- 
ent well widths (see Fig. ^. Notice that the peak of 
the correlation function for the electron-electron couple 
in X~ slowly shifts towards smaller distances as the well 
width decreases (see Fig. ^(a)), at the same time the 
tail of the function becomes smaller. The peak of the 
electron-hole correlation function also increases (see Fig. 
^(b)) but is still centered around zero. With decreas- 
ing well width the X~ becomes less extended. A similar 
behavior was observed for X+. 

In Fig. ^ we show the contour plots of 
\4'oi~^ie,^2e,^h)\'^ for a negatively charged exciton in 
a quantum well of width 100 A, where lighter regions cor- 
respond to lower probability. In Figs. 0(a,b) we plot the 
projection of the electron probability density in the xy- 
and xz-plane when the hole is fixed at T^/i — (0, 0, 0) and 
one of the two electrons is fixed at T*ie = (1.5aB,0,0). 
The distance between the particles is equal to the aver- 
age electron- hole distance in the X^. The two symbols 
show the positions of the two fixed particles. The second 
electron sits close to the hole and the fixed electron sees 
an exciton consisting of the hole-second electron couple. 
Notice, from Fig. 0(b), that the second electron slightly 
penetrates the barrier. In Figs. 0(c,d) we plot the pro- 
jection of the probability density in the xy- and xz-plane 
when the hole is fixed at 'r'h — (0,0,0) and one of the 
two electrons is fixed at T^ie = (as, 0,0). Thus, the dis- 
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tance between the electron and the hole is now smaller 
than the average distance. Now, in the xy-plane a small 
part of the second electron sits on top of the hole and the 
largest part is situated outside the white ellipse defined 
by the position of the fixed electron. Thus now the fixed 
electron and hole act like an exciton to which the second 
electron is bound. The situation is similar in the xz-plane 
where part of the second electron sits on top of the hole 
and the other part is almost symmetrically distributed in 
two puddles around x=±2.2as. 

In Fig. p| wc show the contour plots of 
0o(^i?i, ^2?i, ^e)P for a positively charged exciton in 
a quantum well of width 100 A. In Figs. ^(a,b), we plot 
the projection of the electron probability density func- 
tion on the xy- and the xz-plane where we fixed the two 
holes at r^i;, = (0, 0, 0) and 1^211 = {2.2aB, 0, 0). The dis- 
tance between the holes is now equal to their average dis- 
tance. Notice that the electron is now equally distributed 
over the two holes. Notice also (see Fig. |8|(b)) that the 
electrons do not penetrate into the barrier, as opposite 
to what happens for the X~, indicating that the elec- 
tron is now more strongly bound. In Fig. ||(c), we show 
again the projection of the probability density function 
on the xy-plane, where now the second hole is fixed at 
~r'2h — (1.5a_B, 0, 0). In this case the lobes of the electron 
probability density function repel each other and are no 
longer centered on the position of the holes. The holes are 
closer than their average distance but the distance among 
the centers of the lobes is still approximately equal to the 
hole-hole average distance in the X+. In Fig. ||(d) we fix 
the position of the electron 'r'e — (0, 0, 0) and the first 
hole 'r'lh = (1.5a_B, 0, 0), such that the average distance 
is equal to the electron-hole pair average distance. The 
second hole is now centered around the electron which is 
the same as the behavior of the second electron in X~ (see 
Fig. 0(a)). We found that the contour plot of the proba- 
bility density of finding the hole in a X~ at a position 
when the two electrons are fixed, is practically the same 
as the one for the electron in the X+ with the two holes 
fixed (see Fig. ||(a)). 

In Fig. H we show the wavefunction of the hole in X^ 
and the electron in X+ along the cc-axis when the two 
particles having the same charge are fixed (solid circles in 
Fig. ^). Notice that although the trions are in a singlet 
state the wavefunction is anti-symmetric for reflections 
around the mid-point between the two fixed particles. 
The reason is that interchanging the two fixed particles 
must result in a sign change of the wavefunction. Re- 
mark that the electron in X+ is more localized on the two 
holes as compared to the hole in X~ which is spread out 
over the two electrons. However in both cases the wave- 
function has a node between the two fixed particles, in 
contrast to what would happen if the interaction among 
the particles would be "chemical bonding" -like. It seems 
then reasonable to say that in the same way in which the 
X^ can be described as an exciton with an extra elec- 
tron moving around the electron-hole couple and weakly 
bound to it, the X+ can be viewed as an exciton in which 



an extra hole moves in a orbit around the electron-hole 
couple. The latter picture is different from a system in 
which two holes bind through an electron, i.e. H2'-like. 
Another confirmation for this picture comes from the pair 
correlation functions. Suppose that the electron, for X^, 
is in the origin, then the hole will be near the origin as 
indicated by the electron-hole correlation function and 
the other electron will be situated around the position 
of the peak in the electron-electron pair correlation func- 
tion. So the picture we get is an electron-hole pair with 
an extra electron moving around it. If we switch the role 
of hole and electron, a similar picture can be imagined 
for the X+, with the only difference being that now the 
extra hole sits even further from the electron-hole couple 
than in the X~. Thus, the charged exciton is similar to 
the charged positron. The similarity in the structure of 
the two different species of charged excitons is consistent 
with the fact that their correlation energy is found to be 
equal. 

IV. COMPARISON WITH EXPERIMENTS 

Experimental data in zero magnetic field wese reported 
for the_binding energy OiL-the X~ in a 100 Ata, 200 m , 
220 AB and in a 300 AE3 quantum well. The reported 
values are 2.1, 1.15, l.liO.l and 0.9±0.1 meV respec- 
tively and are compared in Fig. |l^ with our theoretical 
results. The value of = 2 meV for a 80 A well is for 
a GaAa/AlAs quantum well and was measured by Yan 
et aZ.Ell. The theoretical results for the binding energy 
are represented by a shaded region which gives the ac- 
curacy of our calculation for the binding energy. Note 
that the accuracy obtained for the total energy is bet- 
ter than 1%, however, its error propagates and increases 
because of the subtractions (see Eq. (||)) which have to 
be made in order to obtain the binding energy, which is 
one order of magnitude lower than the total energy. An 
important consequence of this observation is that any ap- 
proximation made in the calculation of E{X^) may lead 
to substantial errors in the binding energy. For compar- 
ison we also report (open squares) thp_theoretical result 
obtained by Tsuchiya and KatayamaE3 using the quan- 
tum Monte Carlo method. Notice that the results of 
Ref. ^ agree very well with ours, however our calcula- 
tion goes down to smaller well width. The binding en- 
ergy first increases with increasing quantum well width 
and then, after reaching a maximum of -EB=1.6meV at 
L « 35 A starts to decrease. The decrease becomes very 
slow for quantum well width above 70A. The increase of 
the binding energy with decreasing quantum well width 
agrees qualitatively with the experimental data, but the 
experimental increase is much faster than the one we 
find theoretically. The inclusion of the conduction band 
non-parabolicity would increase the binding energy only 
slightly. We believe that the increased discrepancy be- 
tween theory and experiment with decreasing well width 
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is a consequence of the localization of the trion due to the 
presence of quant uiHrjpejll width fluctuations, as was also 
found for biexcitonal3E3. This is consistent with the fact 
that for L=300A our result agrees with the experiments 
and that the effect of the quantum well width modulation 
on the localization of the exciton and the trion increases 
with decreasing well width. 

A similar calculation was done for CdTe quantum 
well structures (Fig. ^(a)) and ZnSe based quantum 
well structures (Fig. pd[(b)). The binding energy ver- 
sus the well width in these materials is shown in Figs. 
|TT|(a-h)-..(splid curve) and is compared to experimental 
datcO'ESE3 (symbols in Figs. |ll|(a,b)). The parame- 
ters used in the calculation for CdTe based structures 
are Ve = 216.4 meV, Yh = 163 meV, me=0.096mo, 
m^=0.19mo, as = 54 A, which results into Rj^ = 13.8 
meV. The value of the barriers are taken from Ref. |2j. 
Notice that for this structure we have the same potential 
barrier heights than for the GaAs case, however the ratio 
between the mass of the electron and the mass of the hole 
is very different, namely me/m^ = 0.505 (CdTe) as com- 
pared to me/m/i — 0.196 (GaAs). In the range 200A <L 
<600 A the theoretical curve is shifted by about 1 meV 
with respect to the experimental results. Below 200A the 
experimental results increase faster with decreasing L as 
compared to the experimental data which is probably a 
consequence of the above mentioned increased localiza- 
tion of the trion. 

For the ZnSe structure we use the parameters of the 
ZnSe/ZnBcMgSc structures, AV= 230 meV with Ve = 
0.70AV, Yh = 0.30AV, m^ = 0.16mo, m,, = 0.8mo, 
as = 30.05 A, which results in Rj, — 53.34 meV. 
The results are qualitatively similar to the one obtained 
for the GaAs/AlGaAs quantum wells. The agreement 
with experiment is also in this case not satisfactory (see 
Fig.[ri|(b)) except for the 200 A wide quantum well. 

To understand the fact that the theoretical results for 
CdTe- and ZnSe-based structures underestimate so much 
the experimental data even for large well widths, we have 
to take into account that these materials are strongly po- 
lar. In the present work we are neglecting polaronic ef- 
fects and it is known, at least for the case of excitons, 
that this leads to sja underestimation of the binding en- 
ergy of the systemcj. Currently only a calculation of the 
polaron copection to the ground state of a system 
is availablecZI but no calculation for the trion system has 
been published. For the D~ system we know that the po- 
laron correction equals the 3D polaron correction down 
to rather small well widths. When, we shift the results 
in Figs. 0(a,b) by the constant values l.lmeV and 1.6 
meV, respectively (dotted curve in Figs. |ll](a,b), they 
agree very well with the experimental results over a large 
range of quantum well widths. We belipie that these 
shifts are due to polaron effects. Shi et aZO obtained an 
upper limit of « QAahw^o to the polaron contribution 
to the binding energy of the D~ system in a wide quan- 
tum well (a is the electron-phonon coupling constant and 
hjJLO is the optical phonon energy). In an X~ system the 



hole is not localized which will strongly reduce the po- 
laron effect to an estimated value of 0.1-0.2 ahiOLO- For 
CdTe quantum wells with a — 0.3 andJiujLo — 21.1 meV 
this gives 0.6-1.2 meV while for ZnSeES we have a — 0.42 
and hujLo = 31.5 meV and consequently 1.3-2.6 meV. 
These value are comparable to the shifts in Fig. [ll|(a,b), 
and agree with the fact that the shifts for the ZnSe quan- 
tum wells is larger than for CdTe quantum wells. 



V. CONCLUSION 

In this paper we applied the stochastic variational 
method to study the ground state of the exciton and the 
charged exciton in a quantum well. This is the first time, 
to our knowledge, that a calculation fully includes the 
effect of the Coulomb interaction and the confinement 
due to the quantum well, and thus the particle-particle 
correlation in both the direction of the quantum well and 
the confinement direction. The results obtained do not 
show a big qualitative difference from the one already 
present in the literature, however substantial quantita- 
tive differences are found. This difference leads to an 
improvement in the agreement with experimental data. 
However, the experimentally measured binding energy 
for a negatively charged exciton increases much faster 
with increasing well width at small well width than our 
theoretical results. We believe that this discrepancy is a 
consequence of the increased localization of the exciton 
and trion with decreasing well width. A similar|-cpnclu- 
sion was also reached recently for biexcitonsB30. For 
CdTe- and ZnSe-based quantum wells the polaron effect, 
which was not included in our approach, is expected to 
lead to a substantial shift (^ 1 meV) of the binding en- 
ergy to larger values. Also in this case the trapping of the 
trions on the quantum well width fluctuations is proba- 
bly responsible for the rapid increase of the trion binding 
energy below Lr:! 100 A . The study of the conditional 
probability distribution of the particles in the system and 
of the pair correlation functions lead us to conclude that 
a charged exciton is similar to a charged positron. This 
conclusion is important as it supports the fact that the 
correlation energy for and X+ is found to be equal. 
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FIG. 1. The correlation energy of the exciton and the neg- 
ative charged exciton vs. the quantum well width. 



8 



> 
B 

a" 




GaAs/AIGaAs 


G a A S] 


s 






hole / 
















electron 




-60 -40 -20 





20 


40 60 




z (A) 







no mass-mismatch 
mass-mismatch 



20 30 40 50 60 70 
L(A) 



FIG. 2. The correlation energy of the negative charged ex- 
citon vs. the quantum well width, for the case of constant 
masses in the well and in the barrier, and for the the case of 
different electron and hole masses in the well and in the bar- 
rier of a GaAs/Alo.sGao.TAs quantum well. In the inset the 
wave function for both the electron and the hole are shown. 



FIG. 3. The total energy of the negative charged exciton 
and the positive charged exciton vs. mh/rrie for a 200 A wide 
quantum well and for a quantum well of width 100 A (inset) . 
The total energy of a in the same quantum well is given 
by the dash-dotted line for comparison. 




FIG. 4. The correlation energy of the negatively charged 
exciton (symbols) and the one of the positively charged exci- 
ton (sohd curve) vs. the well width. The correlation energy 
of the exciton (dotted curve) is given as reference. 
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FIG. 5. 3D pair correlation function for different pairs of 
particles in and in vs. the distance between the par- 
ticles. In the inset the 2D pair correlation function is shown. 
The curve convention is the same in the two plots. 



FIG. 6. The pair correlation function for the elec- 
tron-electron (a) and electron-hole (b) in X~ for different 
quantum well widths . 




(a): 



I 



(b) 



-0.6 
-0.8 



FIG. 7. Contour maps of the conditional probability for the 
X~ in a quantum well of width lOOAsi as. The fixed particles 
are indicated by symbols (circle with a cross for the hole and 
circle with a minus sign for the electron). The dotted lines 
indicate the quantum well boundaries. 
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FIG. 10. The theoretical (shaded curve) and experimental 
(jjj(open circles) binding energies of the negatively charged exci- 
^on in a GaAs/AlGaAs quantum well vs. the well width. The 
theoretical results of the quantum Monte Carlo calculation of 
Ref.22 are shown by the open squares. 



FIG. 8. Contour maps of the conditional probability for 
the X+ in a quantum well of width 100 Awas. The fixed 
particles are indicated by symbols (circle with a cross for the 
hole and circle with a minus for the electron). The dotted 
lines indicate the quantum well boundaries. 
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FIG. 9. The wavefunction of the hole in when the two 
electrons are fixed (dotted curve) and of the electron in 
when the two holes are fixed (dashed curve) along the di- 
rection [1,0,0]. The two solid dots indicate the two equally 
charged fixed particles. 
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FIG. 11. The binding energy of the negative charged exci- 
ton vs. the quantum well width for the charged exciton for 
CdTe-based structures (a) and for ZnSe-based structures (b) . 
The experimental data axe taken from Refs. 2,24, for the 
CdTe-based structures, and from Ref. 25 for the ZnSe-based 
structures. Our theoretical results are given by the solid 
curve. The dotted curve is our result shifted by a constant. 
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